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Abstract 

We study effects of both magnetic field and radiation on unsteady 
free convective magnetohydrodynamic (MHD) flow under the influ- 
ence of heat and mass transfer. An implicit finite difference scheme 
of Crank-Nicolson method is used to analyze the results. The gov- 
erning boundary layer equations along with the initial and boundary 
conditions are casted into a dimensionless form which results in a 
more precise solution. We determine the change in velocity under the 
influence of radiation and magnetic field with temperature and con- 
centration taken into account. A code is constructed (FORTRAN 6.5) 
to obtain predictions following from the numerical solution which have 
been compared with the existing information in the literature, and a 
very good agreement is obtained. The results obtained through this 
study will be helpful in the prediction of flow and heat transfer in the 
presence of magnetic field and radiation. 
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1 Introduction 



The flow problems occur in many engineering, technology and space sciences 
applications. From an industrial point of view, the radiation effect of free 
convective magnetohydrodynamic (MHD) flow with heat transfer has become 
more important in recent years. Many processes in engineering phenomena 
occur at high temperatures and the knowledge of radiation heat transfer 
has become very significant for the design of pertinent equipment at these 
temperatures. The unsteady natural convection flow over a vertical cylinder 
with MHD flow has been given an extensive attention in literature. Sev- 
eral theoretical and experimental investigations on heat and mass transfer 
characteristics, radiation effects and free convection have been carried out 
from the beginning of sixties. Heat transfer problem from vertical circular 
cyl inders to the surrounding in itially quiescent fluids has first been analyzed 
by iGoldstein and Briggsl (119641 ) for the transient free convection. The effects 
of heat and mass transfer o n the naturally c onvec ti ve flow along a verti- 



cal cy linder were studied by IChen and Yuh I (119801 ) . iRaptis and Kafousias 



(119821 ) have reported the influence of magnetic field upon the steady free 
convective flow through a porous medium bounded by an infinite vertical 
plate with constant section velocity where t he plate temperature w as as- 
sumed constant. By using harmonic analysis, iTakhar and Rami (ll99ll ) have 
focused the effects of Hall currents on MHD free convection boundary layer 
flow via a porous medium past through a plate. IVelusamy and Gars (119921 ) 
have given a numerical solution for the transient natural convection over 
heat by generation of vertical cylinders of various thermal capacities and 
radii. However, a fully implicit finite difference tech n ique w a s used t o solve 
non- linear set of equations by IVelusamy and Gard (1l992f ) . iHelmyl (119981 ) 
studied the unsteady laminar free convection flow of an electrically conduct- 
ing fluid through a porous medium bounded by an infinite vertical plane 
surface of constant temperature. The radiati on effect on h eat tr ansfer over 
a stretching surface h a s first been studied by lElbashbeshy (l2000h . Recently, 



Ghaly and Elbarbaryl ( 120001 ) investigated the radiation effect on MHD free 



convective flow of a gas at a stretching surface with a uniform free stream. 



Seddeekl (120011 ) has studied thermal radiation and buoyancy effects on MHD 



free convective heat generated flow o yer an accelerating perme a ble su rface 
with temperature dependent viscosity. lAboeldahab and Seddeekl (120011 ) have 
studied the radiation effects on unsteady MHD free convection Hall current 
near and infinite vertical porous plate. 
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However, no magnetic field and radiation working together on unsteady 
free convective flow along vertical cylinder under the influence of heat and 
mass transfer has been reported. Due to this fact, the main objective of the 
present work is to investigate the problem of unsteady free convective MHD 
flow of an incompressible viscous fluid passing through a vertical cylinder un- 
der the influence of a magnetic field and radiation. Thus, under the influence 
of radiation and magnetic field we were able to achieve the important changes 
in velocity, temperature and concentration. For this purpose, results from 
an implicit finite difference scheme of Crank-Nicolson method were used. In 
§2.1, the governing equations are given. In §2.2 the details of the numerical 
solutions are described and the main features of the flow are presented. In §3, 
the numerical results and discussions are presented. In §4, the conclusions 
are given. 



2 Numerical database 
2.1 The governing equations 

In this section, the effect of radiation on unsteady free convective MHD 
flow is analyzed. The governing equations are solved more accurately than 
those in the previous works on the subject. The numerical predictions 
have been compared with the existing info r mation in the l i teratu re (e.g. 
Pullepu. Ekambavanan and Chamkha ( 2008 ); Deka and Paul ( 201 ll ) ) and a 



very good agreement is obtained. Axial coordinate x is measured along the 
axis of the cylinder. The radial coordinate r is measured normal to the re- 
axis of the cylinder. We consider the MHD flow of a viscous incompressible 
fluid past through a vertical cylinder of radius tq with transversely applied 
magnetic field. The radiating gas is said to be non-gray if its absorption 
coefficient is dependent on wavelength. The equation that described the 
conservation of radiative transfer in a unit volume for all wavelengths is 



DC 



V • q r = J k x (T) (4vr/ bA (T) - G x ) dX (1) 
o 

where q r is the average heat flux of the radiation, is the spectral intensity 
for a black body, k\ is the absorption coefficient, and the incident radiation 
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G\ is defined as 

g x = f i bx (n)dn (2) 

Q=47T 

where Q is the solid angle. For an optically thin fluid exchanging radiation 
with an isothermal flat plate, according to Equation ([2]) and Kirchhoff law, 
the incident radiation is given by 

G x = 4irI bX (T w )=4e bX (T w ) (3) 

where T is the average temperature value of the porous plate. Equation (JTJ 
is then reduced to 

oo 

V • q r = J k Xw (T)(e bX (T) - e bX (T w )) d\ (4) 
o 

where k Xw = k x (T w ) is the mean absorption coefficient, e bX is Planck's func- 
tion, and T is the temperature of the fluid in a boundary layer. After ex- 
panding e bX {T) and k Xw {T) in Taylor series around T w for a small parameter 
(T — T w ), the Equation (j3J) becomes 

V -q r = -AT(T-T w ) (5) 

where T is a constant of integration. All fluid properties are considered to be 
constant except for the density variation which induces the buoyancy force. 
The viscous dissipation is neglected in the energy equation. The governing 
equations are expressed in terms of conservation of mass, momentum, en- 
ergy, mass diffusion and the interaction of the flow with the magnetic field. 
We assumed that there is no chemical reaction between diffusing species and 
the fluid. Initially, it is assumed that the fluid and the cylinder are of the 
same temperature and also at same concentration C'^. Then the surface 
temperature and concentration of the cylinder are raised to a uniform tem- 
perature and concentration and C' w , respectively. Here, the subscripts w 
and oo indicate the conditions on the wall and free stream, respectively. It 
is further assumed that the interaction of the induced axial magnetic field 
with the flow is considered to be negligible compared to the interaction of 
the applied magnetic field B . Then the flow is governed by the following 
system of equations 
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(6) 

(7) 
(8) 
(9) 



The equations above indicate the continuity, momentum and energy con- 
servation, and mass diffusion equations, respectively. Here, u, v, t', g, B, 
T', B*, C, a, p, a and D represent the x-component of velocity, the radial 
component of velocity, time, acceleration due to gravity, magnetic field (for 
temperature), temperature, modified magnetic field (for concentration), con- 
centration, Stefan-Boltzmann constant, density, thermal diffusivity of fluid 
and mass diffusion coefficient, respectively. The initial and boundary condi- 
tions are given by 



t' < : 
V 



V 

OO) 



f > : 

r 



v 



U 

a 

u 

a 



o, 

C' 



c' w 



v = 0, 

for all x and r 

v = 0, 
at r = r 



(10) 



We introduce the following dimensionless quantities 
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X = Gr 1 x/tq, R = r/r , U = ur Gr 1 /v, 

V — vro/v, t — vt'jr\, Sc = v/D, 

F = AYrl/{vGy/r), N = Gr*/(Gr), M = aB$r$/(pv), 

rpl rpl 

Pr = v/a, a = k/(pc p ), 9 



W OO 

c C'-C'n Cr _ gBr*{r w - T'J ^ _ gB*rl(C' w - C'^ 



C' W -C'J 



(11) 



Here X, R, U, V, t, N, M, Pr, Gr, Gr*, 6, Sc, k, c p and F are the dimen- 
sionless variables such as spatial coordinate, radial coordinate, x-component 
of velocity, radial velocity, time, the combined buoyancy ratio parameter, 
magnetic field parameter, Prandtl number, Grashof number (for tempera- 
ture), modified Grashof number (for concentration), temperature, Schmidt 
number, thermal conductivity, specific heat at constant pressure and the ra- 
diation parameter, respectively. The Equations ©-([9]) are reduced to the 
following form 

86 TT dQ Tr d0 Id f„d6\ „,„ , , A , 

m +u W + v lm = l^m{ R m) +F(e - l) (14) 
sc ,,dc „ac i 8 I ae\ 



The corresponding initial and boundary conditions expressed in the non- 
dimensional quantities are given by: 
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if < : £7 = 0, 1/ = 0, 



t' > : 



= 0, (7 = 0, for all X and i? 



[7 = 0, V = 0, 

= 1, C = l, 

£7 = 0, = 0, 

£7^0, 0^0, 



at R = 1 
C = at X = 
C ->■ as i? -)■ oo 



(16) 



2.2 Numerical technique 

The unsteady, non-linear, coupled partial differential equations f jT2|) - fri5l) 
with the initial and boundary conditions f fl6|) are solved by employing a finite 
difference sch eme of Crank-Nicholson method which was d is cussed by many 
authors (e.g. iPullepu. Ekambavanan and Chamkhal (120081 ); iDeka and Paul 



( 2011 )). The Crank- Nicholson method is a finite difference method (ICrank and Nicolson 
( 119471 )) and it solves the heat equation and similar partial differential equa- 
tions numerically. For the heat equation, the method consider the normalized 
heat equation in one dimension with homogeneous Dirichlet boundary con- 
ditions. This method is based on central difference in space and gives the 
second-order convergence in time. Equivalently, it is the average of forward 
and backward Euler schemes in time. In the applications to large time steps 
or high spatial resolution, the less accurate backward Euler method can be 
used. The Euler me t hod is both stable and independent from oscillations 
( ICrank and Nicolson! ( 119471 )). The partial differential equations (}l2"]) - (jl5p 
are converted into finite difference equations: 



— (£77^ -UZQ + UtfA- £7"^ + £7£ - UUj + £7*-! - £7« 1>i _ 1 ) 
+ 2KR {V ^ 1 ~ V ^~ l + *3 " V ^~ l] = ° (17) 



1 

At 



m 1 -m+ ^- x u u u if ui hj ) 
i 



_i T/n+l/rrn+l _ rrn+1 , rrn _ jjn \ 

AAR ' j ' -.i i ' ' • i u i,i-i> 
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+ 2(KW iU ^ ~ 2U ^ + U ^ + U ^ - 2U ^ + (18) 

^m 1 - <%) + ^tm l - + %j - out) 

-n\m 1 +^- 1 )) (19) 



~ + 2Ax ' j ' j ~ J J ~ 



4A# 



J (-C?/ 1 - Cg& + 2C£ +1 - (20) 



25c(Ai?) 2 



The code is constructed by using FORTRAN 6.5 to obtain the numerical 
solutions of the finite difference equations which were obtained from the 
Crank-Nicolson method. 



3 Results and discussions 

In order to get an insight into the physics of the problem, we present the re- 
sults of numerical computations of velocity, temperature and concentration 
for different values of Prandtl number Pr, Schmidt number Sc, magnetic 
field parameter M, combined buoyancy ratio parameter N, time t and radi- 
ation parameter F. Pr physically relates the relative thickness of hydrody- 
namic and thermal boundary layers while Sc relates the relative thi c kness of 



hydrodynamic and concentration boundary layers (iDeka and Paul! (120111 )). 
Velocity profiles represented by Figfj] show the effect of Pr and Sc numbers 
for different values of time near the cylinder at X = 1. Effects of Pr and Sc 
numbers on the transient temperature profiles are given in Figj2] and tran- 
sient concentration profiles in Figj3] with respect to time, near to cylinder 
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at X = 1. The effect of Pr number is very important in temperature field 
especially. The radiation effect had not been taken into account in the first 
three figures. 

The transient profiles of velocity, temperature and concentration in the 
presence of radiation parameter F are indicated in FigJU Figj5] and Figj6] 
with respect to spatial coordinate X, magnetic field parameter M, combined 
buoyancy ratio parameter N, time t, Prandtl Pr and Schmidt Sc numbers. 
All these numbers are in the dimensionless form as indicated above. To study 
the behavior of the velocity, temperature and concentration profiles, a code 
using FORTRAN (version 6.5) is constructed to solve Equations ([17]) - pD]) . 
Following these results from the simulation for an unsteady free convective 
MHD flow past a vertical cylinder with heat and mass transfer, the following 
general points can be made: 

i) In Figfj], the velocity increases with time and reaches the steady state in 
a certain lapse of time. Temporal maximum is not observed. The R direction 
denotes the radial coordinate and U direction denotes the x-component of the 
velocity. The Schmidt number Sc is taken as 0.7, while the Prandtl numbers 
Pr are indicated as 0.7 and 7.0 with solid and dashed lines, respectively. 
Hence, the transient velocity decreases with increase in Pr number. When 
N, the combined the ratio parameter, increases the velocity increases. The 
velocity profile decreases with increasing magnetic field parameter M. Time 
taken to reach the steady state increases with increasing M. 

ii) In FigJ2l higher Pr values show higher heat transfer rates. The R direc- 
tion denotes the radial coordinate and U direction denotes the x-component 
of the velocity without the effect of radiation. The transient temperature 
profiles increase with decrease in the value of Pr but increases with increas- 
ing time. The higher values of time t show the steady state conditions. The 
Prandtl numbers Pr are indicated as 0.7 and 7.0 with solid and dashed lines, 
respectively. The temperature profiles increase with increasing magnetic field 
parameter M but it will have decreasing slopes with increasing temperatures 
in the limit of the combined buoyancy ratio parameter N —> 1. 

iii) Unsteady concentration profiles are shown in FigfSJ The R direction 
denotes the radial coordinate and U direction denotes the x-component of 
the velocity. We hold the Prandtl number Pr steady and change the Schmidt 
number Sc. When Sc increases, the mass transfer rate increases. The solid 
and dashed lines represent different values of Sc numbers as 0.7 and 7.0, 
respectively. The transient concentration decreases with the increase in value 
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of Sc and increases with increasing time t and reached the steady state for 
larger time t. The concentration profiles increase with increasing magnetic 
field parameter M while it does not show larger effect in the limit of the 
combined buoyancy ratio parameter N — >■ 1. In the case of increasing N, the 
concentration profiles are observed to have decreasing slopes with increasing 
concentration. 

iv) The effect of radiation parameter F on the velocity, temperature and 
concentration is shown in FigJU Figj5] and Figj6j respectively. It is observed 
that the velocity, temperature and concentration increase as the radiation 
parameter F increases. This result qualitatively agrees with expectations 
since the effects of radiation and surface temperature are to increase the 
rates of energy transport to the fluid, thus increasing the temperature of the 
fluid. 

v) In FigJU the variation of velocity profiles in the presence of radiation 
parameter are indicated. The R direction denotes the radial coordinate and 
U direction denotes the x-component of the velocity without (F = 0) and 
with (F = 2, 4, 5, 6) the effect of radiation. The velocity profiles (and slopes) 
increase with increasing radiation values. The velocity is related to the rate 
of energy transport of the fluid. As long as the radiation increases the energy 
transport rate, the velocity will also be increased. 

vi) There is a direct relationship between temperature and radiation. It 
can clearly be seen in FigJSJ The R direction denotes the radial coordinate 
and T direction denotes the temperature without (F = 0) and with (F = 
2, 4, 5, 6) the effect of radiation. Therefore, the temperature of the fluid 
increases with increasing radiation and increasing the radiation will stop 
the sudden decline in the temperature as seen in the slope of the bottom 
temperature profile. 

vii) The effect of radiation for the concentration profiles is seen in FigJHJ 
The R direction denotes the radial coordinate and C direction denotes the 
concentration without (F = 0) and with (F = 2, 4, 6) the effect of radiation. 
Although there was not a very big difference in concentration, the noticeable 
increase in the radial coordinate is very obvious. Mostly, the radiation effects 
the velocity and temperature profiles. The less deep but similar profiles can 
also be seen in concentration. 
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1 2 R 3 4 5 



Figure 1: The velocity profiles for M and iV at different times with the effects 
of Pr numbers, Pr = 0.7,7.0 and Sc number, Sc = 0.7. Transient velocity 
profiles are indicated at X = 1 and F = 0. The value of i with star (*) symbol 
denotes the time taken to reach steady state. The R direction denotes the 
radial coordinate and U direction denotes x-component of the velocity. Here, 
M and iV indicate the magnetic field parameter and the combined buoyancy 
ratio parameter, respectively. 
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Figure 2: The temperature profiles for M and iV at different times with the 
effects of Pr numbers, Pr = 0.7,7.0 and Sc number, Sc = 0.7. Transient 
temperature profiles are indicated at X = 1 and F = 0. The value of t with 
star (*) symbol denotes the time taken to reach steady state. The R direction 
denotes the radial coordinate and T direction denotes the temperature. 
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Figure 3: The concentration profiles for M and N at different times with the 
effects of Pr numbers, Pr = 0.7, 7.0 and Sc number, Sc = 0.7. Transient 
concentration profiles at X = 1 and F — 0. The value of t with star (*) 
symbol denotes the time taken to reach steady state. The R direction denotes 
the radial coordinate and C direction denotes concentration. 
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C. E 




Figure 4: Transient velocity profiles at X = 1, M = 4, N = 1, t = 1, 
Sc = 0.7 and Pr = 0.7. Different plots show different radiation effects for 
the velocity profiles. From the bottom profile to the top, the values of 0, 2, 4, 5 
and 6 indicate the changes of temperature without and with radiation, re- 
spectively. Increasing radiation values have increased the slopes of velocity 
profiles. The R direction denotes the radial coordinate and U direction de- 
notes x-component of the velocity. 



14 




H 1 1 1 =H 

1 2^3 4 



Figure 5: Transient temperature profiles at X — 1, M = 4, N = 1, t = 1, 
Sc = 0.7 and Pr = 0.7. Different plots show different radiation effects for 
the temperature profiles. From the bottom profile to the top, the values of 
0,2,4,5 and 6 indicate the changes of temperature without and with radia- 
tion, respectively. In the absence of radiation, a very deep slope stands out. 
The R direction denotes the radial coordinate and T direction denotes the 
temperature. 
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Figure 6: Transient concentration profiles at X — 1, M — 4, N — 1, t — 1, 
Sc = 0.7 and Pr = 0.7. Different plots show different radiation effects for 
the concentration profiles. From the bottom profile to the top, the values of 
0, 2, 4 and 6 indicate the changes of concentration without and with radiation, 
respectively. The R direction denotes the radial coordinate and C direction 
denotes the concentration. 



16 



4 Conclusions 



In this work, we have made a numerical analysis for unsteady free convective 
MHD flow past a vertical cylinder with heat and mass transfer under the 
effect of magnetic field and radiation by employing finite difference scheme 
of Crack-Nicolson method. For the first time, we included in our analysis the 
magnetic field and radiation parameters together. The velocity, temperature 
and concentration profiles are drawn for different values of parameters such 
as Pr, Sc, M, N, t and F to study the behavior of velocity, temperature and 
concentration of the flow starting from the initial radial coordinate P = 4. 
Present results have been compared with the available results in literature 
and are found to be in a good agreement. Later, we have examined the 
effect of radiation for the same profiles. The radiation will be more effective 
on the profiles of velocity, temperature and concentration. Increasing the 
radiation implies increasing of velocity and temperature in the flow. The flow 
problems at high operating temperatures where the radiation effects are quite 
significant occur in many engineering applications, science and technology. 
The effects of free convective heat flow play an important role in agriculture, 
petroleum, industries, geological formations. The results obtained from this 
study will be helpful in the prediction of flow and heat transfer. 
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